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Abstract 

Angular distributions of ct(6'), Ay, iTu, T20, T21, and T22 have been mea- 
sured for d-p scattering at Ec.m. = 667 keV. This set of high-precision data is 
compared to variational calculations with the nucleon-nucleon potential alone 
and also to calculations including a three-nucleon (3N) potential. Agreement 
with cross-section and tensor analyzing power data is excellent when a 3N 
potential is used. However, a comparison between the vector analyzing pow- 
ers reveals differences of approximately 40% in the maxima of the angular 
distributions which is larger than reported at higher energies for both p-d 
and n-d scattering. Single-energy phase-shift analyses were performed on this 
data set and a similar data set at Ec.m. = 431.3 keV. The role of the different 
phase-shift parameters in fitting these data is discussed. 
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I. INTRODUCTION 



An important question in low-energy nuclear physics is whether a nucleon-nucleon (NN) 
interaction is sufficient to describe few-nucleon systems or whether a three-nucleon force 
(3NF) is necessary when the system consists of A> 3. To solve for the bound and scattering 
states of the three-nucleon (3N) system, state-of-the-art theoretical methods like the Fad- 
deev method (which works well for n-d scattering) have been developed. It was not until 
the advancement of variational techniques which could include the Coulomb force rigor- 
ously that serious comparisons could be made with p-d scattering data. Current theoretical 
investigations incorporate phenomenological NN potentials such as CD Bonn 0, Nijmegen 
III Q, and Argonne vis (AV18) 0. These potentials are based on the one-pion-exchange 
potential for large interparticle separations as well as fits to the present NN database. When 
calculations using those models are performed in the three-nucleon system, the ^He and 
binding energies are underestimated by approximately 0.5 to 1.0 MeV depending on which 
NN potential is employed. The under-prediction of the A = 3 binding energies is the first 
indication that the two-nucleon interaction alone, as parameterized by the NN potentials in 
their current forms, is insufficient for an accurate description of the 3N data. To remedy 
the situation, phenomenological 3N potentials such as the Tuscon- Melbourne (TM) 0, the 
Brazil (BR) and the Urbana IX (UR-IX) |^ were introduced which are based on the 27r 
exchange interaction and adjusted to reproduce the triton binding energies. By applying 
chiral-perturbation theory to the 3NFs, Friar et al. showed that the present phenomeno- 
logical 3NFs are equivalent to ffist-order in the Lagrangian they constructed except for one 
term in the TM potential which was found to be of next-to-leading order. 

Beyond the binding-energy problem, there is another large discrepancy in the properties 
of the 3N system, namely the vector analyzing powers (VAPs) Ay and iTn for N-d scattering. 
Comparisons between variational calculations |ffOl including one of the NN potentials and 



VAP measurements at Ecm. < 2 MeV |llT|JT^ show that the calculations under-predict the 



data by ~ 30% in the maximum of the angular distributions. The disagreement becomes 



worse at lower energies where the difference increases to ~ 40% at Ecm. = 431.3 keV |13,14 
This VAP problem was ffist observed for Ay in n-d scattering and has been labelled the 
''Ay puzzle" [p!5| , p!6| . The trend of the ''Ay puzzle" for both Ay and iTn is a decreasing 
discrepancy as the energy increases as shown in Fig. |T|. The disagreement mostly disappears 



sX Ecm ~ 30 MeV 17 . The inclusion of a 3NF in the variational calculations lessens the 



discrepancies with Ay and iTu data for p-d scattering at the aforementioned energies by 
no more than 15%. The effect of the magnetic moments was investigated by Stoks [|1^] and 
found not to provide the solution. 

To investigate this question further, we have measured with high precision Ay and iTn 
as well as the cross section a{0) and the tensor analyzing powers (TAPs) T20, T21, T22 
for p-d elastic scattering at Ecm. = 667 keV. To isolate the deficient phase shifts and 
mixing parameters in the theoretical calculations, two single-energy phase-shift analyses are 
performed, one with these data at Ecm. = 667 keV and another with the same observables at 
Ecm. = 431.3 keV previously measured by our group. Additional information 



concerning the experimental details and analysis is available in Ref. [^. Moreover, the 
667 keV cross-section measurement was included in a study with other cross-section 
data below Erm = 2 MeV. 



2 



II. EXPERIMENTAL DETAILS 



The basic experimental setup and techniques were developed for our measurements of d-p 

= 431.3 keV and are described in Ref. lT|Jl|Jl|-|| 



elastic scattering at Ec 



The modifica- 
= 667±1 keV 



tions to the setup described in Ref. |jT^ for the experiments conducted at Ecm 
will be discussed. One advantage of making measurements at Ecm. = 667 keV over 431.3 keV 
is the reduction of the Rutherford scattering of the incident particles from carbon which was 
present in the targets. The statistical uncertainties in this set of measurements are signif- 
icantly smaller than in the measurement at 431.3 keV. The only slight disadvantage is the 
presence of additional proton groups from the ^^C((i, p)^^C reaction in the deuteron beam 
experiments. All six angular distribution measurements employed the FN tandem accelera- 
tor at the Triangle Universities Nuclear Laboratory to give a deuteron (proton) beam energy 
in the center of the target of 2.00 MeV (1.00 MeV). 



A. Cross-Section Measurements 

The relative cross-section was measured with a beam of deuterons on thin hydrogenated 
carbon targets ^3] which contained approximately (0.5 ± 0.1) x 10^*^ H/cm^ and (1.0 ± 
0.2) X 10^^ C/cm^. The rest of the experimental details can be found in Ref. A 
sample spectrum for ^H(c?, (i)^H scattering at 6iab = 26.1° is shown in Fig. |^. The absolute 
cross-section was obtained by normalizing the relative a (6) for d-p scattering to the well- 
known ^H(p, p)^H cross-section determined by the Nijmegen phase-shift analysis |2^. The 



uncertainties of the absolute cross-section measurement are 0.5% from counting statistics, 
0.5% from the normalization to the p-p scattering cross-section, and 0.6% from the beam- 
current integration. The absolute a (6) angular distributions are shown in Fig. The data 
are divided by the results of the variational calculation to emphasize the differences between 
the data and the calculations. 



B. Analyzing Power Measurements 



For the five analyzing power measurements, the techniques used were similar to those 
described in Refs. ^j^. The targets [||] contained (1.0 ± 0.2) x 10^^ hydro gen or deu- 
terium atoms and (2.0 ± 0.4) x 10^^ carbon atoms per cm^. The beam polarization for each 
experiment was measured online with a polarimeter situated behind the scattering chamber. 
For the T20, T21, and T22 experiments, a polarimeter ||2^ based on the ^He((i,p)^He reaction 
was utilized. The beam polarizations for the T20 and T22 measurements were pzz ~ ±0.8 
and Pzz ~ for the three spin states with a 3% error. For the T21 measurement, the beam 
polarizations were Pzz ~ ±0.8 and Pzz ~ with a 2% error. For iTu, a polarimeter ||21 



based on the ^H((i, p)^H reaction was used, and the beam polarizations were Pz ~ ±0.55. 
For the proton beam, the ^Li(p, ^He)^He reaction pT] , |26| was employed for polarimetry. The 
beam polarization over the duration of the experiment was pz ~ ±0.65. The error in the 
beam polarizations for both experiments was 2%. 

The coincidence technique described in Ref. [|I^ was used for the VAP experiments. A 
new addition was 2-/im-thick mylar foils that were placed in front of the detectors to block 
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any heavy recoil particles such as carbon atoms. By making the measurements at the higher 
energy, it was possible to cover a larger angular range {Oc.m. ~ 50° — 140°) than obtained at 
431.3 keV. A typical timing spectrum for p-d coincidences for a deuteron beam incident on 
a hydrogen target is shown in Fig. ^. The TAP and VAP angular distributions are given in 
Figs. I and ||, respectively. 



III. COMPARISONS WITH THE THEORETICAL CALCULATIONS 



Theoretical calculations to compare with the data have been produced using the AV18 
potential and the AV18 plus a 3NF. These calculations employ the Pair-Correlated Hyper- 
spherical Harmonic basis to construct the scattering wave functions and the Kohn variational 
principle to solve for the scattering-matrix elements . In order to study the sensitivity to 
different 3NF models, the UR-IX, the TM and its modified version denoted TM' interactions 
have been considered. Originally the TM potential ^ was given as a sum of four terms char- 
acterized by the strengths a, b, c, d. To assure consistency with chiral symmetry the a-term 
in the TM model has been redefined with a new constant a' = a — 2m\c (which essentially 
implies a change of sign for the constant a) and the c-term has been set to zero |^ . The TM' 
model, as well as the TM and the UR-IX potentials, have been used recently in an extensive 
study of 3NF effects in n-d scattering above Ecm. = 2 MeV ^\ . For completeness, here we 



will present results using the same set of constants as Ref. |23] for the TM and TM' models. 

Accompanying the data presented in Figs. |^ and ^ are several separate variational cal- 
culations, one including the AV18 potential and the others with the AVIS and either the 
UR-IX or TM potentials. It is clear that the calculations with or without the 3N poten- 
tial reproduce (y{6) and the TAP data fairly well. On the other hand, the agreement with 
the VAP data is much poorer. These results are qualitatively similar to those observed 
at higher energies which were mentioned in the Introduction. For the present case, the 
differences between the NN-potential-only calculations and the Ay and iTn measurements 
are approximately 40% and 37%, respectively at the maximum of the angular distributions. 
With the inclusion of the UR-IX potential, the agreement for all of the observables improves; 
however, the change in the VAPs is only marginal and a large discrepancy still remains on 
the order of 36% and 29% for Ay and iTu, respectively. When the TM 3N-potential is 
included, the VAP calculations show slightly worse agreement with the data (39% and 34% 
differences) than in the case of the UR-IX potential. However, when the TM' potential is 
used, the calculations are closer to the VAP data than for the case of the UR-IX potential; 
the difference between Ay{iTii) data and the calculations is 32%(25%). The same trend has 
been observed in Ref. at Ec.m. = 2 MeV. The difference between the calculations with 
the TM' potential and those with UR-IX is about 5%. While each of the 3NFs mentioned 
here provide some improvement to the Ay and iTn calculations, none of the potentials elim- 
inates the descrepancy. To obtain a more quantitative comparison, the reduced x^, X%i was 
calculated for each observable, and the results are described in Sec. |V. 

Besides the TM' potential, other modified NN and 3N force models have recently become 
available. For example, a new NN interaction constructed from chiral perturbation theory 
seems to give a better description of the VAP data at low energy ||2^. A different possibility 



may be in the construction of a new 3N force. Hiiber and Friar |29] have suggested that a 



possible candidate should be a spin-dependent 3NF since the VAP is a difference between 
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polarized cross-sections. Following this direction, a phenomenological 3N L ■ S force has 
recently been proposed [^]. This new force essentially modifies the scalar function v[\{r) 
already present in the L ■ S term of the AV18 potential in triplet spin and isospin channels. 
Explicitly the following form has been proposed 

= E^n(^u)L., ■ S,,Pn(u) + W^oe-"''ELy- ■ S,,Pn(^J) • (1) 

i<j i<j 

where Pu is a projector in channels with spin Sij = 1 and isospin Tij = 1. The hyperradius 
p is 



p = y 3(^12+^23+^31) (2) 

and Wq and a are parameters characterizing the strength and range of the three-body term. 
This force influences the splitting of the ^Pj phases and the magnitude of €3^2- without 
appreciably affecting the other parameters. Moreover, when the deuteron and spectator 
nucleon are separated by a large distance, the proposed interaction reduces to the original 
NN L ■ S potential v[\{rij). Fig. ^ shows Ay and iTu calculated using this new potential 
for three different strength and range parameters. The values for Wq and a have been 



taken from Ref. and have been fixed in order to reproduce the VAP data at Ec.m. = 2 
MeV. The addition of the 3N spin-orbit force has provided much better agreement with our 
VAP data than previous calculations. The various phase-shift parameters will be discussed 
in detail in Sec. [V|. Moreover the calculated differential cross sections and TAPs remain 
essentially unchanged by the new term. 



IV. CALCULATIONS OF REDUCED 

To investigate quantitatively the effects of the UR-IX potential, reduced calculations 
were completed with both the AV18 potential and the AV18+UR-IX potentials. More details 
for the a{6) comparison are presented in Ref. and the results are listed in the first two 



lines of Table |. None of the data were renormalized in order to find a minimal Xn- From xff 
calculations, it is possible to draw conclusions similar to those from the visual comparisons: 
the a{6) and the TAP data are reproduced very well while VAP data show poorer agreement 
with theory. However, all of the x'n improved when a 3NF is included. For a{6), x'n 
improves by an order of magnitude. The x'n is reduced significantly for each TAP. The 
calculations of Ay and iTn with UR-IX show improved agreement with the data; however, 
the x'n remains in the hundreds. 

V. PHASE-SHIFT ANALYSIS 

The technique of Phase Shift Analysis (PSA) is useful for determining the depen- 
dences of the observables on the various partial waves, and also for studying the sources 
of disagreements between theory and experiment. Here we consider data at two energies: 
Ec.m. = 667 keV, the energy of the present data, and also Ec.m. = 431.3 keV, the energy of 



previous measurements in our laboratory [jl4[. The PSAs reported here are carried out in 
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the framework of Seyler , where the scattering matrix is parameterized in terms of phase 
shifts and mixing parameters. The phase shift parameters in general depend on orbital 
angular momentum L, total angular momentum J, and the channel spin S. The mixing 
parameters connect partial waves with different L and/or S for a given J and parity. 

Before embarking upon a detailed PSA, it is important to understand the convergence 
of the observables versus the angular momentum cutoff. Studies were performed using the 
AV18+UR-IX calculations for the phase shifts and mixing parameters. The dependence on 
the angular momentum cutoff L^ax was examined by setting to zero all phase shifts and 
mixing parameters with L > L^ax for different values of L^ax- Note that the effect of 
the Coulomb potential is retained in all partial waves, so it is expected that the observ- 
ables will converge quickly. We have quantified the convergence of the observables by using 
Xn comparisons of the calculations with different Lmax values to the experimental data at 
Ec.m. = 431.3 and 667 keV. The results are shown in Table First of all we find very little 
change in x'n between L^ax = 3 and L^ax = 4, indicating, as expected, that the results 
are already well-converged with L^ax = 3 at these low energies. It is also interesting to 
note that the cross section and TAPs are fairly well described with Lmax = 2, while only 
Lmax = 1 is required to produce a similar level of convergence for the VAPs. Smaller, but 
statistically-significant improvements are obtained by extending to Lmax = 3. We have also 
observed in these studies that the observables T20 and T22 are highly dependent on the L = 2 
phase shifts and mixing parameters. 

Single-energy PSAs have been performed by using the phase shifts and mixing parameters 
produced by the AV18-t-UR-IX potential as starting values. Some of the parameters are 
allowed to vary freely; the best-fit values for the variable parameters are determined by 
minimization with attention given to the fitting procedure to avoid local minima. Phase 
shifts were included up to a maximum angular momentum of L = 4, but only those with 
L < 2 were allowed to vary. 

Previous phase-shift analyses pi]|JTT| ] have linked the YAP problem to the P-wave phase- 
shifts. To test the effects of the different P-wave phase-shifts and their associated mixing 
parameters, two-parameter fits where made. Tables Hand III show the x'n ^wo energies 

for sample trials where only the P-wave phase-shifts were varied. By varying the ^Pj phase- 
shifts and/or the mixing parameter 61/2" freely (trials 1 and 2), there is very little change 
in the x'n with the exception of iTu at 431.3 keV. With the combination of £3/2- aiid a ^Pj 
phase-shift (trials 3, 4, and 5), good fits were produced with the most dramatic improvement 
occurring in each VAP. At E^m. = 667 keV, the x'n reduced by a factor of approximately 
300 for Ay and a factor of 30 for iTu. For the E^m. = 431.3 keV case, the fits produced 
reductions in x'n by factors of approximately 20 and 10 for Ay and iTu, respectively, with 
the exception of the fit to iTu with the free parameters ^P3/2 and 63/2-- A comparison of 
the three ^Pj trials indicates that the two-parameter fit with the ^Pi/2 phase-shift and £3/2- 
gives the best result although the differences with the other two trials are very small. These 
trials so far illustrate the influence of the P-wave phase-shifts and 63/2-! however, a more 
realistic fit (trial 6) is to allow all of the S- and P-wave phase-shifts with e^/2- ^^^1 63^2^ 
vary freely (9 parameters in all). The remaining phases and mixing parameters were taken 
from the AV18-I-UR-IX calculations [jlO[- The results of the fits are listed in Table 

While there were improvements with the full S- and P-wave fits over the two parameter 
^Pj fits, the reduction in x'n small. To make a more quantitative comparison between 



6 



trials, errors were determined for specific parameters by the change in a specific parameter 
necessary to increase the minimal by 1 [^. Allowing more parameters to vary led to 
a larger error in each parameter which underlies our attempts to find the best fit with the 
least number of parameters. The uncertainties in phase shifts and mixing parameters when 
all S- and P-wave parameters were varied are given in Table 0. These errors reflect the 
contribution of statistical uncertainties in the data only. 

The PSA presented here produced very good fits to the data sets at Ec.m. = 431.3 keV and 
667 keV with sensible results. The greatest influence on the Ay and iTn fits came from the 
^Pj phase-shifts and £3/2- parameter. This result was noticed in the earlier analyses pT!|JTT| 
at higher energies. It is difficult to provide a more precise value from the PSA for each phase- 
shift with the present data sets at 431.3 keV and 667 keV. What these analyses do confirm 
is the need to modify the ^Pj phase-shifts and the mixing between the ^P3/2 and the ^P3/2 
phase-shifts (as evident from the change in 63/2-)- In fact, €3/2- was increased (in absolute 
value) by about 19% from the 667 keV PSA and by about 28% from the 431.3 keV PSA 
over the predicted value even when the UR-IX potential was included in the calculations. 
The reader should be cautioned that the errors from both analyses allow much wider ranges 
of values for £3/2- • The theoretical calculations also appear to underestimate the splitting 
AP = '^P5/2— For example at 667 keV the calculated splitting is AP = 0.87° compared 
to 0.96° obtained in the fit. At 431.3 keV the difference is even bigger; the calculated value 
is AP = 0.48° compared to 0.63° from the fit. The small, but nevertheless important, 



underestimation of AP has been found also in the calculations ai Erm =2 MeV 10 



VI. CONCLUSIONS 

In conclusion, we have presented high-precision measurements of <j{9), Ay, iTn, T20, T21, 
and T22 for p-d elastic scattering at a center-of-mass energy of 667 keV. The cross section 
and TAP data show excellent agreement with the variational calculations using the AV18 
and UR-IX potentials. The "Ay puzzle" persists at this energy with the inclusion of the UR- 
IX in the calculations giving a small improvement in the agreement with the VAP data. In 
order to extend the analysis to different forms of the 3NFs, the AV18-I-TM and AV18+TM' 
potential models have also been considered in the comparisons to the data. The calculations 
with either TM potential under-predict the VAP data by a similar amount as in the case of 
the AV18+UR-IX results. From the reduced analysis, it is clear that adding a 3NF into 
the calculations improves the agreement with the data for all of the observables. 

From a single-energy PSA at very low energies (£'c.m. = 431.3 keV and 667 keV), the 
under-prediction of the VAP data by the calculations appears to be caused by certain P- 
wave components. By varying two parameters, 63^2- and one of the '^Pj phase-shifts, and 
maintaining all of the other parameters as determined from theory, an almost perfect fit 
was achieved for Ay. For iTn, the Xn reduced by more than an order of magnitude to 
approximately 1.3 for the 431.3 keV PSA and 1.8 for the 667 keV analysis. To obtain an even 
better fit, all of the S- and P-wave phase-shifts along with e^/2- 63/2- were varied freely 
to obtain Xn of ~ 1.0 for each observable. The results of this analysis are the same within 
the errors as the PSA with 63^2- ^'^d a ^Pj phase-shift as the only free parameters. From 
comparisons of the PSA results to calculations it is evident that the theory underestimates 
the splitting of the ^Pj phase shifts as well as the mixing between ^P3/2 and ^P3/2- 
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The inclusion of a phenomenological 3N spin-orbit force, which fixes the P-wave mixing 
in the calculations, provides better agreement with the VAP data without disturbing the 
agreement obtained for the cross-section and TAP data. While this fact does not determine 
the existence of a 3NF, it indicates a need to develop and refine the current 3N potential 
models p8[. In addition, one can not rule out the possibility that currently employed NN 



28 



Pj interactions are incorrect 

Finally, the Ay problem can be studied with information from p-^He elastic scattering 
data. The variational calculations employed in this paper have recently been extended to the 



p- He scattering observables pPI . In this case there is an underprediction by approximately 



40% of new Ay data taken at an -Ec.m. below 2 MeV. Moreover, Ay for p-^He scattering is 
an order of magnitude larger than for the p-d scattering case. This fact makes it an ideal 
candidate for future investigations, and more p-'^He elastic scattering data at low energies 
is needed. 
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TABLES 



TABLE I. The Xn calculated for each trial of the phase-shift analysis for each observable 
measured at Ec.m. = 667 keV. The number of data points in the angular distribution of each 
observable is the value in parentheses under the title. The total Xn sum of for each 

observable divided by 143 (the total number of data points). 





Free 








Xn 








Trial 


Parameters 


a{e) 


Ay 


iTii 


T20 


T21 




Total 






(56) 


(7) 


(8) 


(24) 


(24) 


(24) 


(143) 


AV18 


None 


45.2 


275.8 


112.4 


3.5 


3.5 


2.8 


39.3 


AV18+UR 


None 


1.2 


190.5 


61.4 


1.0 


2.5 


0.7 


13.9 


1 


-'^1/2) -^3/2 


1.1 


199.8 


53.3 


0.9 


1.9 


0.7 


13.8 


2 


ei/2- 


1.2 


187.7 


46.3 


1.9 


4.9 


0.8 


13.5 


3 


'^Pl/2, £3/2- 


1.1 


0.6 


1.8 


1.7 


1.2 


0.7 


1.2 


4 


^-^3/2) £3/2- 


1.2 


0.5 


5.3 


2.2 


5.7 


0.7 


2.2 


5 


'^^5/2) £3/2- 


1.3 


0.5 


1.9 


0.8 


1.8 


0.7 


1.2 


6 


^'S'l/2) '^Pj, 

'^Pj^ ei/2-) £3/2- 


1.0 


0.8 


1.1 


0.8 


0.9 


0.5 


0.9 



TABLE IL The convergence of the partial wave expansion is indicated by x'n fo'^ different 
values of Lj^ax at Ecm. = 431.3 and 667 keV. The number of data points for each observable is 
shown in parentheses. 





a{e) 


Ay 




x% 

T20 


T21 


^22 


Total 


Ecm. = 431.3 keV 




(22) 


(3) 


(3) 


(24) 


(22) 


(20) 


(94) 




1 


8.45 


17.97 


14.16 


42.15 


18.08 


29.73 


24.32 




2 


2.31 


17.58 


13.17 


1.53 


5.40 


1.72 


3.54 




3 


1.87 


17.62 


13.13 


1.91 


2.61 


1.55 


2.85 




4 


1.86 


17.62 


13.12 


1.94 


2.61 


1.53 


2.85 


Ec.m. = 667 keV 




(56) 


(7) 


(8) 


(24) 


(24) 


(24) 


(143) 




1 


49.22 


194.82 


69.74 


118.49 


51.34 


202.61 


95.22 




2 


1.37 


189.92 


61.88 


2.44 


10.98 


3.96 


16.21 




3 


1.11 


190.43 


61.38 


1.11 


2.83 


0.66 


13.96 




4 


1.15 


190.46 


61.39 


0.98 


2.52 


0.66 


13.91 
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TABLE III. The Xn calculated for specific trials of the phase-shift analysis at 
Ec.m. = 431.3 keV. The number of data points in the angular distribution of each observable 
is the value in parentheses under the title. The total Xn is the sum of for each observable 
divided by 94 (the total number of data points). 





Free 








9 

Xn 








Trial 


Parameters 


aie) 


Ay 


iTu 


T20 


721 




Total 






(22) 


(3) 


(3) 


(24) 


(22) 


(20) 


(94) 


AV18 


None 


23.2 


23.5 


19.3 


2.9 


4.2 


2.1 


9.0 


AV18+UR 


None 


1.9 


17.6 


13.1 


1.9 


2.6 


1.5 


2.8 


1 


2 p 2 p 

^1 /2 ) ^3/2 


1.6 


25.1 


2.7 


1.3 


1.1 


1.6 


2.2 


2 




1.6 


30.1 


2.2 


1.2 


0.9 


1.6 


2.2 


3 


f3/2 


1.7 


0.9 


1.3 


1.2 


0.9 


1.5 


1.3 


4 


''-^3/2) ^3/2- 


1.7 


0.5 


16.6 


2.0 


2.5 


1.5 


2.4 


5 


4 p 

-^5/2) ^3/2- 


2.1 


1.4 


3.4 


1.6 


1.9 


1.5 


1.8 


6 


^5i/2, '^'S'3/2) ^Pj, 
^Pj> ^1/2-: ^3/2- 


1.1 


0.3 


0.9 


1.1 


0.8 


1.3 


1.0 



TABLE IV. Phase shifts and mixing parameters obtained from single-energy phase-shift anal- 
ysis where all S- and P-wave phase-shifts were varied. The first column for each energy contains 
the values calculated with the AV18 and UR-IX potentials and were the starting values for the 
analyses. The second column lists the values produced from the fits to the data. 





E — 


431.3 keV 


Pern. 


= 667 keV 


Paramol cr 


C'alculali()n(") 


Fii n 


C'ak;ulalion(°) 


Fil C^) 


^'S'1/2 


-5.811 


-7.64 ±0.74 


-10.00 


-10.53 ±0.59 


'^'S'3/2 


-28.08 


-27.48 ±0.25 


-36.93 


-36.81 ±0.19 


2 p 

-^1/2 


-1.993 


-2.17 ±0.44 


-3.339 


-2.91 ±0.15 


2p 

-^3/2 


-1.975 


-1.61 ±0.28 


-3.303 


-3.42 ± 0.13 


^1/2- 


2.499 


1.82 ± 0.46 


2.954 


3.33 ±0.13 


4p 


5.262 


5.06 ±0.11 


9.199 


9.18 ±0.07 


'P3/2 


6.127 


6.19 ±0.12 


10.69 


10.62 ±0.05 


'P^I2 


5.743 


5.69 ±0.10 


10.07 


10.14 ±0.04 


€3/2- 


-0.886 


-1.13 ±0.29 


-1.060 


-1.25 ±0.06 
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FIGURES 
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FIG. 1. Differences between experimental data and theoretical calculations using the AV18 
potential at the maximum of the angular distributions for (a) iTn and (b) Ay. The points are 
plotted as the difference between experiment and theory divided by the experimental value. The 
data are taken from the following sources: Ref. 14j20[| (circles), Ref. 
Ref. |11,35| (open squares), Ref. (triangles), and Ref. ||3^ (diamond). 
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FIG. 2. Typical spectrum for a 2.00 MeV deuteron beam on thin hydrogenated carbon foil. 
The detectors were positioned at Oiab = 26.1°. In addition to the peak at channel 360, there are 
two smaller peaks at channels 160 and 210 which arise from ^^C(d,p)^^C reactions. 
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FIG. 3. Angular distributions for p-d elastic scattering at Ec.m. = 667 keV. Graph (a) shows 
the ratio of the cross-sections to the AVIS+UR-IX calculation. The circles represent the data, 
including statistical errors, divided by the AV18+UR-IX calculation. The solid line is the ratio of 
the AV18 calculation to the AV18+UR-IX calculation. Graphs (b), (c), and (d) show the angular 
distributions of the TAP data. The solid and dashed curves are calculations with the AVIS and 
AV18+UR-IX potentials, respectively. The dot-dashed line marks zero on these graphs. 
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FIG. 4. A proton-ductcron time-of-flight spectrum using a four-detector configuration. The 
time resolution is approximately 20 ns. 
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FIG. 5. Angular distributions for iTn and Ay for p-d elastic scattering at Ecm. = 667 keV. 
The errors include the uncertainty in the beam polarization as well as statistical uncertainties. The 
solid, dashed, and dotted curves are variational calculations with the AV18, AV18+UR-IX, and 
AV18+TM potentials, respectively. The dot-dashed curve is the calculation with the AVIS+TM' 
potentials. 
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FIG. 6. Comparison between iTn and Ay data and calculations with AV18+UR-IX with the 
addition of phenomenological spin-orbit forces at Ec.m. = 667 keV. The solid curve is the calculation 
without any modifications. The dotted curve represents a long-range interaction (Wq = — 1 MeV, 
a = 0.7 fm~^), the dashed curve is for a medium-range {Wo = -10 MeV, a = 1.2 fm^^), and the 
dot-dashed curve indicates a short-range interaction {Wq = -20 MeV, a = 1.5 fm"^) 
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